The Search For A Permanent Electric Dipole Moment Using Atomic Indium 



B. K. Sahoo * 

Theoretical Physics Division, Physical Research Laboratory, Ahmedahad-380009, India 



o 

(N 

o 
O 

(N 



R. Pandey and B. P. Das 
Theoretical Astrophysics Group, Indian Institute of Astrophysics, Bangalore-560034, India 
(Dated: Received date; Accepted date) 

We propose indium (In) as a possible candidate for observing the permanent electric dipole mo- 
ment (EDM) arising from the violations of parity (P) and time-reversal (T) symmetries. This atom 
has been laser cooled and therefore the measurement of its EDM has the potential of improving on 
the current best EDM limit for a paramagnetic atom which comes from thallium. We report the 
results of our calculations of the EDM enhancement factor due to the electron EDM and the ratio 
of the atomic EDM to the electron-nucleus scalar-pseudoscalar (S-PS) interaction coupling constant 
in In in the framework of the relativistic coupled cluster theory. It might be possible to get new 
limits for the electron EDM and the S-PS CP violating coupling constant by combining the results 
of our calculations with the measured value of the EDM of In when it is available. These limits 
could have important implications for the standard model (SM) of particle physics. 
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It is now widely recognized that atomic electric dipole 
moments (EDMs) arising from the violations of parity 
(P) and time-reversal (T) symmetries can provide impor- 
tant information about new physics beyond the standard 
model (SM) [H, T violation implies CP violation via 
the CPT theorem 3. The dominant sources of the EDM 
of a paramagnetic atom are the EDM of an electron and 
the scalar-pseudoscalar (S-PS) interaction between the 
electron and the nucleus which violates P as well as T 
symmetries Q. Atomic EDMs due to the electron EDM 
and the S-PS electron-nucleus interaction can shed light 
on CP violation in the leptonic and semi-leptonic sectors 
[U, Q . The origin of both of these kinds of CP violations 
are not well understood. 

The best limit on the EDM of a paramagnetic atom 
currently comes from thallium (Tl) [6- 9]. A new gener- 
ation of EDM experiments on the alkali atoms like ru- 
bidium [l3|, caesium [lol-|T2| and francium [l^ based on 
the techniques of laser cooling and trapping are currently 
underway. These experiments in principle have the ad- 
vantages of both the beam and the cell experiments [l3| . 
In particular it is possible to apply large electric fields 
and the coherence times are long in these experiments 
p^ . The projected precision of the current alkali EDM 
experiments are about two orders of magnitude better 
than that of the Tl experiment [lol - [l^ . Indium (In), 
which is homologous to Tl appears to be an attractive 
candidate for the search of a permanent EDM for essen- 
tially two reasons: First, this atom has been laser cooled 
[isj and second, the EDM enhancement factor due to the 
electron EDM and the ratio of the atomic EDM to the S- 
PS coupling constant can be calculated more accurately 



than that of Tl due to its relatively smaller size. To 
demonstrate the second point more elaborately, we have 
carried out ab initio calculations for the corresponding 
factors due to the EDM of electron and the S-PS electron- 
nucleus interaction in In and discuss the role of different 
correlation effects in these properties here. 

Following the work of Sandars [3] and its extension 
Tie -18], the effective atomic EDM Hamiltonian due to 
the electron EDM can be written as 
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and the S-PS interaction Hamiltonian is given by 
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S-PS 
EDM 
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where de is the intrinsic e-EDM, 7^ is a pseudo-scalar 
Dirac matrix, C5 is the dimensionless S-PS constant and 
Pn{fj) is the jth electron density over the nucleus. 

The above interaction Hamiltonians mix atomic states 
of opposite parities but with the same angular momen- 
tum. They can be treated as first order perturbations 
as their strengths are sufficiently weak. Therefore any 
atomic state with a valence electron 'v' can be expressed 
after the inclusion of these interactions as 



|vI/„) = |vI'(o))+A|*W), 
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where \^v) is the modified wave function to the original 
wave function j^i"'') by the first order correction l^'i^''). 
A represents the weak coupling parameter for H'^j^j^j 

or Cs for H|Bm ■ 

The enhancement factor due to the electron EDM and 
the ratio of the atomic EDM to the S-PS coupling con- 
stant which is denoted hy TZ = ^ here for the atomic 
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EDM Da of a state is given by 



7^ = 
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We consider only single and double excitations in the 
expansion of the RCC wave functions (CCSD approxi- 
mation), by defining, 



T — T1+T2 and 5„ = 5*11, + S2V , 



(0.10) 



where D is the electric dipole (El) operator. With the ^^^j^ perturbed and unperturbed operators. Fur- 

explicit form of j'l'l, ), we get 



ther, we construct triple excitation operators for 5'^°^ as. 
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where / represents the intermediate states and Hedm is 
one of the EDM interactions given above. 

The above expression depends explicitly on El matrix 
elements, excitation energies (EEs) and matrix elements 
of Hedm- However, it is possible to consider only a fi- 
nite number of excited states in the evaluation of this 
quantity if a sum-over-states approach is used. In or- 
der to circumvent this problem, we solve the first order 
perturbed wave function in a similar approach of 



(ij(")-£;("))|vi/W) = -if^^M|*r>- 
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In the above expression, ij'"^ is the atomic Hamiltonian 



and is the energy for the state l^'l,"''). It is possible 
to estimate the accuracy of the corresponding TZ values 
by calculating the properties required to determine Eq. 
(|0.5p for the dominant intermediate states. 

We employ the coupled-cluster theory in the relativis- 
tic framework (RCC theory) to evaluate |vE'i°'') and I^P^"^') 
as described in d, H, [13, [11] . These wave functions can 
be expressed as 



(0)\ 



and 



rp(0) 



rp(0 



{1-f 



(0.7) 

l + 5l°M +4'^}|$.)(0.8) 



where 1$^,) is the Dirac-Fock (DF) wave function ob- 
tained by appending the valence electron v to the closed- 
shell ([4(i^"] 5s^ in the present case) reference wave func- 
tion, r(0) and Si°^ are the excitation operators for the 
core and valence electrons in the unperturbed case, where 
as, T^^'^ and si^^ are their first order corrections. The 
atomic wave functions are calculated using the Dirac- 
Coulomb (DC) Hamiltonian given by, 

^0 = ^{cai-pj-|-(/3j-l)TOjc2-|-y„(r,)} 

i 

+ J2Vcinj), (0.9) 

i<j 

where a and /3 are Dirac matrices, p is the momentum 
operator, Vn{r) is the nuclear potential and Vc{r) is the 
Coulomb potential. 
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which are used to evaluate the CCSD amplitudes itera- 
tively. This is referred to as CCSD(T) approximation. 
Here, is the single particle energy of an orbital i. 
Hence, TZ in the RCC theory is given by 



7^ 
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where the dressed operator D^'^) — e^^' ' D e^^ \ D = er 
is the El operator due to the applied electric field and cc 
represents the complex conjugate terms. The procedure 
for the calculation of the above expression is discussed 
elsewhere @,i,[i3[il]. 



TABLE I: Enhancement factor due to electron EDM and the 
ratio of the atomic EDM to the S-PS coupling constant de- 
noted by TZ to the ground state of In due to de (dimensionless) 
and S-PS (in ^£-Cs/j4; A is the atomic mass number) inter- 
actions obtained using DF and RCC methods. 



Source 


DF 


CCSD CCSD(T) 


de 


-49.53 


-82.35 -82.37 


S-PS 


-31.10 


-52.59 -52.60 



The single particle orbitals in our calculations are a 
linear combination of Gaussian type of orbitals (GTOs). 
They are optimized by comparing the energies and the ra- 
dial integrals of these orbitals with those obtained numer- 
ically from GRASP2 [l9|. We have allowed excitations 
from all the occupied orbitals to unoccupied bound and 
continuum orbitals with a maximum energy of ISOOau. 
This space is sufficiently large for the convergence of the 
results of our calculations. Orbitals up to / = 4 were in- 
cluded in the active space after observing that the contri- 
butions from the orbitals with higher angular momenta 
were very small; inclusion of these orbitals would have 
been computationally expensive with little or no effect 
on the overall results. 

In Table HI we present the TZ values for both and 
S-PS interactions at the DF, CCSD and CCSD(T) lev- 
els. There is a significant difference between the DF and 
CCSD(T) results highlighting the importance of strong 
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TABLE II: Contributions to TZ (with same unit as in Table 
U]) from different CCSD(T) terms. D^*^^ represents effective 
one-body terms from DC), terms containing bare D operator 
are the effective two-body terms from £)("', "Others" and 
"Norm" give contributions from other non-hnear terms and 
normahzation of the wave function, respectively. 



Term from de from S-PS 

From DF 

{DHedm)c + cc -17.56 -11.03 

[dHedm)v + cc -31.98 -20.07 

Important RCC terms 

I^It^^'' + cc -39.79 -25.33 

D^°l S^^^ + cc -44.59 -28.24 

^il-^^t'+cc 10.10 6.39 



Nonlinear RCC terms 



T^(l)''' r)(0) e(0) 


-4.12 


-2.62 




3.55 


2.27 


sT'd^1s^IJ+cc 


-2.06 


-1.29 


t[^^^DT^'>^ +CC 


-4.74 


-3.02 




-1.24 


-0.78 




0.86 


0.56 




-0.97 


-0.61 


Others 


0.21 


0.17 


Norm 


0.83 


0.52 



correlation effects for the reported TZ factors of this sys- 
tem. Small difference between the CCSD and CCSD(T) 
results suggests that the contributions from the triple 
excitations are small. We give below contributions from 
various CCSD(T) terms to understand the roles of dif- 
ferent correlation effects in this property. 

In Table HIl we present the contributions from the core 
and the virtual orbitals to the factors TZ at the DF level 
and individual contributions from different CCSD(T) 
terms. The CCSD(T) contributions in this table have 
been classified as important terms referring to terms 
whose contributions are large (Goldstone diagrams for 
these terms are shown in Fig. [T|) and nonlinear terms 
whose contributions are relatively small. Fig. [IJa) rep- 
resenting D^^^t[^^ involves the lowest order DF contri- 
butions due to the core orbitals (Fig. [TJi)) and some of 
the higher order core-polarization correlation diagrams 
(Fig. [ljii,iii) and more). By comparing the DF con- 
tributions from core orbitals and it is obvious 
that these core-polarization correlation effects contribute 
significantly; larger than the DF core contribution. The 
largest contributions come from the D^^'iS^^ term and 
as shown in Fig. (lUb)), it contains the lowest order 
DF contribution due to the virtual orbitals. It also in- 




FIG. 1: Break down of important perturbed CCSD(T) di- 
agrams into lower order many-body perturbation diagrams. 
Lines with single arrow showing up and down wards repre- 
sent the occupied and unoccupied orbitals, respectively. Line 
with a double arrow represent the valence orbital. 



eludes many important core-polarization (Fig. [Uv,vi)) 
and pair (Fig. [Ijvii)) correlation diagrams. However, 
the net correlation contribution due to this diagram is 
not as large as it is from the higher order diagrams in 
£)(o)j'^(i)^ The other important term that has been con- 
sidered is D^'^^s'^^ and it accounts for mainly a particular 
class of core polarization effects. 

We now consider some of the quantitative aspects of 
the correlation effects in the above electron EDM en- 
hancement factor and the ratio of the atomic EDM to 
the S-PS coupling constant. These effects in In are sig- 
nificantly larger than those in the alkalies, but their can- 
cellations are not as severe as in the case of Tl. The 
correlation effects for both the above quantities are one 
and a half times the total DF results. The core correla- 
tion effects are very strong, because of the large overlap 
between the wave functions of the valence 5pi/2 and the 
outermost core 5s orbitals and the small energy differ- 
ence between them. At the DF level, the contributions 
from the virtual orbitals are larger than those of the core 
orbitals. The total contributions from both £i(")t|^' and 

£'(°)S'JJ,'' are comparable. The contribution from 52^-* 
is significant, but with opposite sign. These contribu- 
tions are from the singly excited perturbed states; mainly 
from the [4c?-'^°]5s5p^ state. There are also some signif- 
icant contributions from the higher order RCC terms; 

especially t[''^^ DT!f\ ^j^^^D^^f \ 5^°^ ^j^^ and 

t[^^ 'iMs!^^ through some of the core-polarization and 
pair-correlation effects. Contributions due to the normal- 
ization of the wave function (norm) and terms containing 
non-linear terms in T^^^ and T^^^ are small. 

It is possible to get a sense of the accuracies of the 
individual quantities that appear in Eq. (jO.Sp by com- 
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TABLE III: Comparison of the excitation energies (in cm~^), 
El matrix element and hyperfine structure constants of low- 
lying states in In with the available experimental results and 
all order SD calculations. Uncertainties from our calculations 
are given in the parentheses. 

This Work All order SD [25] Experiment 



Transition 



Excitation energies 



6s 
7s ■ 
8s 
9 s 



5pi/2 
5pi/2 
5pi/2 
5pi/2 



24290(80) 
36217(90) 
40552(100) 
42640(115) 



23747 
35808 
40126 
42238 



24372.956 [20] 
36301.84 [20] 
40637.0 [20] 
42719.0 [20] 



Transition 



El matrix elements 



6s 


-> 5pi/2 


1.91(1) 


1.91 


1.92(8) [21J 


7s 


5pi/2 


0.56(2) 


0.54 




8s 


5pi/2 


0.31(2) 


0.09 




9s 


^ 5pi/2 


0.19(2) 








State 




Ahyp of "^In 






5Pl/2 


2256(30) 


2306 


2282(40) [22] 




6s 


1611(50) 


1812 


1687.2(6) [23] 




7s 


516(30) 


544.5 


541.1(3) [23] 




8s 


234(20) 


240.8 






9s 


106(10) 


128.1 





paring their calculated values with their corresponding 
experimental results. Experimental values for the EEs 
are available up to very high accuracy [1^ and experi- 
mental values of the El matrix elements can be extracted 
from the lifetime measurements of the available s states. 
The matrix elements of Hedm cannot be measured di- 
rectly. However, the accuracy of this quantity can be 
indirectly estimated from the square root of the prod- 
uct of the magnetic dipole hyperfine structure constants 
(.^hyp) of the appropriate states [1, |^. 

We present the results for the EEs, El matrix elements 
and the hyperfine structure constants of the ground and 
excited s states in Table Hill that contribute significantly 
to TZ. Experimental values for the EEs are available for 
most of the excited states, but the hyperfine structure 
constants are experimentally known only for the 5pi/2, 
6s and 7s states. We obtain the El matrix element of the 
6s — >■ 5pi/2 transition as 1.92(8)au by combining the ex- 
perimental values of the lifetime of the 6s state (7.5(7)s) 
[21| and considering the branching ratio as 2:3 of the 
6s — >■ hpi/2 and 6s — > 6^3/2 transitions [l^] (Ref. [1^ can 



be referred for detailed discussion of these results). For 
EEs, the largest uncertainty comes from the Breit inter- 
action followed by the neglected basis orbitals and triple 
excitations. However, most of the uncertainties to the 
El matrix elements and AhypS come from the latter. All 
our results are in good agreement with the experimental 
results. Our calculated results have been compared with 
those obtained by the all order SD method [25], the RCC 
method containing only the linear terms of our CCSD ap- 
proach, in the above table. Our EEs results are in better 
agreement with accurate experimental data than those 
using the all order SD method for all the excited states 
relevant for our present work. The corresponding DF re- 
sults are given elsewhere (see Ref. [ll]) and it is found in 
[25} that the Breit interaction contributes very little to 
this property in this atom. Therefore, the discrepancies 
between the results reported in [25| and the present work 
could be due to the non-linear terms of the CCSD(T) 
method. Our Ahyp for the 6s state, the most important 
excited state in the calculations of 7?., is also more ac- 
curate than the all order SD result. Our results for the 
other quantities agree reasonably well with available ex- 
perimental data. The uncertainties in various quantities 
are estimated by considering the differences between the 
results of the CCSD(T) and CCSD methods as the upper 
limits to the contributions due to the triple excitations, 
neglected relativistic effects; particularly the Breit inter- 
action and omitted higher angular momentum symmetry 
orbitals in the present calculations. After considering all 
possible uncertainties, the enhancement factor due to the 
electron EDM and the ratio of the atomic EDM to the 
S-PS coupling constant for In are estimated to be —82(5) 
and —53(3), respectively. These results are almost five 
and one and a half times smaller than those in Tl [9] and 
Cs lis! respectively, but three times larger than those in 
Rb [18|. Their accuracies can be further improved by 
using the general RCC theory (27| . 

In conclusion, we propose In as a suitable candidate for 
the search of a permanent EDM. Our theoretical stud- 
ies show that accurate calculations of the electron EDM 
enhancement factor and the ratio of the atomic EDM 
to the S-PS coupling constant of this atom are possible. 
The limits for the electron EDM and the S-PS coupling 
constants that can be extracted by combining these fac- 
tors with the measured value of the EDM of this atom 
when it is available, could provide important information 
about the validity of the SM of particle physics. 

We thank H. S. Nataraj for useful discussions. These 
calculations were carried out using the HPC 3TFL0P 
cluster at PRL and the CDAC ParamPadma TeraFlop 
supercomputer. Bangalore. 
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